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Abstract 

We investigate the possibility of introducing the temperature to the one 
dimensional Sznajd model and propose a natural extension of the original 
model by including other types of interactions. We characterise different 
kinds of equilibria into which the extended system can evolve. We deter- 
mine the consequences of fulfilling the detailed balance condition and we 
prove that in some cases it is equivalent to microscopic reversibility. We 
propose a simple definition of the temperature-like quantity that measures 
the size of fluctuations in the system at equilibrium. The complete list of 
zero-temperature degenerated cases is provided. 
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1. Introduction 

The Sznajd model was originally introduced in [l| in order to describe 
the mechanisms of opinion change in the society. The basic element of this 
approach is the "social validation" , a social phenomenon relying on the fact 
that usually two people sharing common opinion have much bigger influence 
on other people in a group than separate individuals. The model itself and its 
modifications have found numerous applications in sociophysics (see jij for a 
recent review), marketing js], 0], finance \5] and politics p, 0, isj. Apart from 
interest in social science, finance and politics, it also brings some new ideas 
to physics, as noted by Slanina and Lavicka |9|. Some theoretical aspects of 
the model were investigated in |9|,[l0l,[ll|. 

The Sznajd model provides a new scheme of interacting between the par- 
ticles. The basic component of the Sznajd model is a lattice of Ising spins. 
Each spin is in one of two states: "up" or "down" , like in the standard Ising 
model. Here the orientation of a given spin designates the opinion of an 
individual described by this spin, e.g. voting for or against in some politi- 
cal context. In the Sznajd model (and in other outflow models in general) 
we have the following dynamics. Within the common framework of random 
updating some chosen spins (their number depends on the variation of the 
model) influence their outer neighbors, as opposed to e.g. Glauber dynam- 
ics, where two spins affect the spin between them. In the first case we speak 
about outflow dynamics, as information flows outwards, and in the second 
case, we have inflow dynamics. Both kinds of models provide microscopic 



description of macroscopic changes in the system, e.g. phase transitions [12 



Some dynamical aspects of one- and two-dimensional versions of the Sznajd 



model were considered in [13|, [ij, ll5| . Interesting scaling behavior of relax- 



ation times for two-dimensional Sznajd model were reported in [16|. A simple 



underlying dynamical framework that drives global dynamics of the system 



was proposed for both Sznajd and Glauber zero-temperature models in [17 

In almost all studies of outflow-type models the rules that govern dy- 
namics are deterministic, and therefore these models are considered as zero- 
temperature. In the social context this means that all people behave in the 
same, predictable way, which is extremely rare in real situations. In gen- 



eral, different types of social response are possible [18[, such as conformity, 
anticonformity, independence and congruence. These terms are related to 
various ways of changing mind in different conditions (e.g. conformity is the 
core of ,, social validation", as a group tries to make an individual following 



2 



the group's opinion). Certainly human behavior is not predictable in 100% 
and with some probability one can expect any kind of social response. There- 
fore it would be reasonable and interesting to seek a model that incorporates 
some diversity or randomness in human activity. On the physical grounds 
this relates to a question: How to incorporate noise, or temperature, to the 
original model in order to allow for randomness inherent in many real sys- 
tems? Such a development of the Sznajd model will allow us to formulate 
some general remarks about types of equilibrium that can be reached by the 
outflow models. 

In the original Sznajd model there is no room for the temperature, and 
so an extension of the model is needed. Here (in the section two) we propose 
such extended model. Our framework comes from taking into account all 
possible configurations of small vicinity of each spin (or individual). It is 



19, 20 



so general, that it contains some previously examined models [13 
as special cases. In order to introduce the notion of temperature in the 
model one should examine the equilibrium of the system first (the standard 
temperature is well defined in equilibrium only). The discussion of that 
issue is provided in the section three. In the next section we propose a 
natural candidate for the temperature in the extended model. Conclusions 
are drawn in the last (fifth) section. In the first appendix we present a 
complete list of special (degenerated) cases of the extended model, for which 
the dynamics always ceases after some finite time. In the second appendix 
we show the detailed arguments leading from the detailed balance conditions 
to the constraints IQ and ffTOl) . 

It should be also noticed that the results presented in this paper are rigor- 
ous and belong to very few contributions that analyze the outflow dynamics 
analytically |9|. 



2. The extended model 

In the original Sznajd model [H of outflow dynamics one considers a 
system of L Ising spins Si = ±1 on one dimensional lattice with periodic 
boundary conditions. In each update two nearest neighbors are selected at 
random. Let us assume that spins i and i+1 were selected. If they are parallel 
{Si = Si+i), their outer neighbors' spins follow the spins of the selected pair 
(that is Si-i = Si+2 = Si). In the case of antiparallel alignment {Si ^ 
S'j+i) two options are considered: either the outer spins acquire antiparallel 
orientation with respect to their closer neighbor {Si-\ = Si+i and Si+2 = Si) 
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or nothing is changed. The latter option was considered to be more natural 
and was often used in the literature, but some other rules (for antiparallel 



alignment) were also used, see e.g. 13|, |l9|, |20|. The dynamics defined in 



this way is undoubtedly zero-temperature, as no random noise is present in 
the system. All initial states approach one of the two final ferromagnetic 
states (all spins up or all spins down) in a finite time. The distribution of 



the relaxation times of this process was investigated in [15 . 

The original model dynamic rules act only on configurations of the type 
. . . y4[i?i?]A . . . (from now on the selected pair of spins, that affect their outer 
neighbours, we will put in square brackets), where A = ±1 and B = =f1 
stands for different spin orientations. In the cases of antiparallel alignment 
of the selected spins and in the case of all four spins parallel in a row nothing 
happens. The natural extension of the model accounts for allowing spin 
changes in these cases. Even though there are some other ways of introducing 
the noise to the Sznajd model (e.g. [^H), the one proposed here seems the 
most natural. 

We consider a ID lattice with periodic boundary conditions (as in the 
original model). There are L Ising spins on the lattice, and the dynamics of 
the model is given by special rules of changing spins in a single update step. 
At each step two consecutive spins are chosen at random, and they influence 
their direct neighbors according to the rules. Let us consider all possible 
configurations of 4 consecutive spins (two middle spins in brackets will control 
the outcome of the update step). We do not consider the external field, so 
our extended evolution rules must be symmetric with the respect to spin 
reversal (+<-)■—). We consider the action of a selected pair independently in 
each direction, so only three spins are to be taken into account (two selected 
and one to be updated). Thus all different possible elementary cases make 
up the following list: ([AAjA, [AAjE, [AB]A and [AB]B). To determine the 
dynamics we have to introduce a vector of probabilities p = (pi,P2,P3,P4) of 
flipping the third spin in the configurations below: 



Pi 

P2 
P3 



[AA]A [AA]B, (1) 

[AA] B ^ [AA]A, (2) 

[AB] A [AB]B, (3) 
: [AB]B -> [AB]A. (4) 



The values of these four constants (or parameters) are crucial in the model 
and govern the dynamics. For convenience we write the complementary 
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probability: p'i = ^ — Pi for i = 1, ...,4, which describes the chance of not 
flipping the spin in a given configuration. 

From the social point of view the parameters pi denotes the intensity 
of specific social behavior: non-conformity {pi), conformity (^2) and more 
subtle cases of changing one's opinion in view of non-conforming nearest and 
next nearest neighbors (ps and ^4). 



3. Equilibrium 

The notion of equilibrium in statistical physics is fundamental. Many 
properties of the systems are defined only in an equilibrium state (or at 
least in local equilibrium, where the changes of macroscopic quantities are 
relatively slow). In many cases however, systems cannot reach equilibrium 
due to e.g. existence of external forces. But the meaning of equilibrium 
is sometimes not clear - there are several non-equivalent definitions of an 
equilibrium in the literature. Let us consider the dynamical system with a 
finite microstate space X. Let qi{t) be the probability that the system at 
time t is in the microstate i El. The most basic definition of an equilibrium 



22| states that the distribution of probability among all possible microstates 



does not vary with time: 

qi{t) = for all i E X. (5) 

This condition can be rewritten using the notion of transition probabilities 
between the states. Let Wnn' be the probability of transition between the mi- 
crostates: n' — )■ n in a single time step (or probability density for continuous 
time). Then the condition states that the probability of leaving the state 
n balances the probability of arriving at the state n from all other states: 



There is also another, more strict approach (see, e.g. ^^), that de- 
fines equilibrium of the system in more rigorous way - it is required that all 
transitions between each two microstates must balance independently. This 
constraint is known as detailed balance condition (or DBC) and reads: 

Wnn'qn' = Wn'nQn (7) 
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for each pair of the microstates n and n' . The direct consequence of ([7]) is 
the absence of net probabihty flow along cyclic transitions among any three 
states i — )■ j — )■ — )■ i (what, in general, can take place in the "basic" 
equilibrium satisfying only the full balance conditions (|6])). 

There is a fundamental interplay between DBC and reversibility. In the 
theory of Markov processes DBC is equivalent to reversibility (but not 
microscopic reversibility, see below). Let us consider a stationary Markov 
process (all finite-dimensional distributions depend on time differences and 
do not change with time translations). A process is said to be reversible 
if its finite-dimensional distributions do not change with the time reversion 
(the transformation t — r — t for some r) . Then the necessary and sufficient 
condition for a stationary Markov process to be reversible is just fullfiling the 
DBC for a collection of positive numbers constituting the equilibrium distri- 
bution of the process. In statistical mechanics the detailed balance appears 
on a basic level. It is known [25] that the second law of thermodynamics 
(nonnegative entropy production in the isolated system) together with the 
detailed balance is equivalent to the condition of microscopic reversibility: 

Wnn' = Wn'n (8) 

stating that the transition probabilities between two states do not depend on 
the direction of the transition. In deriving Boltzmann H-theorem for diluted 
gases DBC appears as an assumption of cross section invariance with respect 



to time inversion. It is also known [22| that DBC is always satisfied for all 



classical physical isolated systems, for which hamiltonian and macroscopic 
observables are even functions of all momenta. 

Even though there is a distinct difference between the microscopic re- 
versibility condition and the detailed balance condition (the first implies the 
second, but not the opposite, in general) many authors use them interchange- 
ably. 

Let us now come back to the extended model and investigate the equilib- 
rium types the system can reach. Since we deal with the finite size system, 
the equilibrium (in the "basic" meaning ^) is always reached in the infinite 



time limit [22[. But there are two general possibilities, depending on the 
choice of the parameters pi, . . . ,p4 of the model. First class is characterized 
by the existence of absorbing states - the dynamics comes to an end after a 
finite time of reaching a microstate that system cannot leave. These cases 
we call degenerated, as the dynamics stops after a finite time. In this situ- 
ation the set Vt of final microstates may consist of a single element, or for 
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other parameters' value may contain many final microstates, but there are 
no allowed transitions between them. The full list of degenerated cases with 
the appropriate values of pi, . . . ,P4 is provided in the Appendix A. 

In all other cases that make up the second class the equilibrium state 
approached by the system has never vanishing fluctuations and the transition 
probabilities between the microstates realising given equilibrium state do not 
vanish. 

It is an interesting question what kind of equilibrium is reached by the 
system depending on the model parameters - or - what kind of reversibility is 
fulfilled by the dynamics of the model. For the equilibrium from the first class 
all transitions between microstates realising given equilibrium state vanish, 
so DBC is satisfied trivially (both sides of the equation ([7]) equal zero). For 
such cases the dynamics is fully reversible at equilibrium in such sense that no 
change can be present - the system is frozen forever in both time directions. 

On the other hand the issue of satisfying DBC in the second class is by 
no means evident. One of the aims of this work is to clarify this matter. 
The details of the reasoning are provided in the appendix B, here we only 
sketch the main idea and gotten results. Since DBC forbids the net flow 
of probability along any cycle of transitions (e.g. i — )■ j — )■ /c — )■ i against 
2 — )■ A; — )■ j — )■ z for any microstates i,j, A; G X), the skillful choice of triples 
k) and their analysis leads to some constraints on parameters pi, . . . ,p4. 
It is argued in the appendix B, that DBC implies the following constraints: 

Pi = P2, (9) 

and 

P3 = Pa- (10) 

In order to remain in the second class we must assume further that at least 
one pair of probabilities does not vanish: 

pi > or p3 > (11) 

(otherwise we have the degenerated case of number from the appendix A). 
In turn assuming ([9]) and (fTOj) leads directly to DBC, since the obtained 
conditions ^ and fIlOp are just the statement of microscopic reversibility. 
Thus the extended model is truly time-reversible only for parameters' val- 
ues satisfying fl9|)- f|TT|) and then time reversibility (or DBC) is equivalent to 
microscopic reversibility. 
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4. The temperature 

After investigating the nature of equilibrium in the extended model we 
now can safely proceed to discuss the notion of temperature in this context. 
As a temperatTirc is usually interrelated to the size of fluctuations in an 
equilibrium state, it would be natural to assign the null temperature to the 
models from the first class with degenerated dynamics, where no changes are 
present at the equilibrium state. 

In statistical mechanics the inverse of the temperature is usually defined 
as an energy derivative of the entropy, but there is no hamiltonian (energy) 
in the model. We can define however a temperature-like parameter T in such 
a simple way that it would measure the mean rate of microstate change at 
equilibrium. To make it system size independent let us normalize it to one: 
we define T as the average number of elementary changes of the microstate 
in a single Monte Carlo step divided by the size of the system. Thus T = 
indicates no changes at all and T = 1 means that the microstate is changing 
at every possible instant of time. 

In general it is difficult to write the way our temperature T depends on 
all pi, . . . ,p4 in a closed form. For all cases of the degenerated dynamics we 
have obviously T = 0. It appears (from Monte Carlo simulations) that in the 
cases realising the stronger equilibrium (these from the second class) with 
Pi = P2 and = Pi the dependence of the temperature on the parameters 
is particularly simple: 

rj. ^ Pi +P2+P3+P4 ^^2) 

The highest temperature T = 1 possible in the model (the highest level of 
noise) is obtained only for the case Pi = P2 = Ps = Pi = ^- We believe that 
the proposed way of introducing the temperature to the outfiow models will 
prove useful and enable to investigate such systems from a new perspective. 

5. Conclusions 

In this paper we proposed an extension of the Sznajd model that allows 
for wider variety of interactions and room for random fiuctuations that are 
observed in real systems. In this extended model we analysed in details 
the nature of possible equilibria reached by the system and discussed in 
this context the question of reversibility of the dynamics. It appeared that 
there are special cases of the model's parameters, for which the dynamics 
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is fully reversible. For the proper description of fluctuations in the model 
we suggested the temperature-like measure that vanishes for all degenerated 
cases and acquires the particularly simple form in the case of the reversible 
dynamics. Applications of such temperature needs further study and we 
believe that this paper will serve as a good starting point in investigating 
outflow models from a new perspective. 

Appendix A. The degenerated cases 

Detailed investigation of the rules leads one to cases, when the dy- 

namics of the system stops after arriving at some absorbing microstate, from 
which the system cannot escape. Such cases we call degenerated, since no 
further nontrivial dynamics is possible. 

Here we present the complete list of degenerated cases (the set VL consists 
of all final microstates that cannot be escaped from): 

0. pi = 0, P2 = 0, P3 = 0, P4 = 

trivial case of no dynamics at all, each configuration is final 

1. pi >0, P2 = 0, ps = 0, P4 = 

Q consists of configurations with all clusters of size one and two, 
e.g. {... + + - + -- + --...} 

2. pi = 0,P2> 0, P3 = 0, P4 = 

the case corresponding to the standard Sznajd model, for which p2 = 1. 
n = { ... + + + + + + + +..., 

...- + - + - + - +...} 

3. pi =0, p2 = 0, ps > 0, p4 = 

Q consists of configurations with all clusters of size bigger than one, 
e.g. {... + + + + + -- + + +...} 

4. pi = 0, P2 = 0, p3 = 0, P4 > 

n = { ...- + - + - + -+...} 

5. pi >0,p2 = 0, p3 > 0, P4 = 

n = { ... + + -- + + --...} 
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6. pi >0,p2 = 0, p3 = 0, P4 > 

n = { ...- + - + - + -+...} 

7. pi =0,P2> 0, p3 > 0, P4 = 

n = { ... + + + + + + + +..., 
... ...} 

8. pi =0,P2> 0, p3 = 0, P4 > 

n = { ... + + + + + + + +..., 



...- + - + - + - +...} 

9. pi =0,P2> 0, P3 > 0, P4 > 

n = { ... + + + + + + + +..., 
... ...} 

10. Pi > 0, P2 > 0, p3 = 0, P4 > 

n = { ...- + - + - + -+...} 

When the lattice size L is an odd number: a) the cases 4, 5, 6 and 10 are 
no longer degenerated, b) in the cases 2 and 8 the antiferromagnetic state 
does not belong to Q. If L is an even number not divisible by 4 the case 5 is 
no longer degenerated as well. 

Let us make a note on the meaning of the above cases on sociological 
grounds. The case no. is not interesting - nothing happens. The next four 
cases (1.-4.) describe situation, when only one type of transition is allowed 
- one kind of behavior: conformist (the case 2. - the Sznajd model), pure 
nonconformist (the case 1.) or more subtle nonconformist (the cases 3. and 
4.). The cases 5.-8. allow for two different types of behavior that compete 
with each other, but always one of them finally wins. The last two cases 
(9. and 10.) admit three kinds of behavior at a time, but finally one of 
them always wins. Apart from the cases 0., 1., 3. and 5. the system ends 
in a stable state of either perfect conformity or total disagreement (each two 
nearest neighbors do not agree with each other). 

Appendix B. Derivation of constraints ([9]) and (1101) 

We now show how satisfying DBC in the extended Sznajd model leads to 
the constraints iQ and ffTOj) . The main idea is to require equality between 



10 



the rates of cyclic transitions between suitably chosen three microstates in 
both directions. 

Let us consider transitions between the following three microstates: 



(/): ... (+) + +(+) ... (B.l) 

(//): ... (+) + +(-) ... (B.2) 

(///): ... (-) + +(-) ... (B.3) 

The dots . . . above describe parts of the configurations that do not contribute 
to the changes and are of no interest for us. 

In the following we calculate transition rates for cyclic changes /—;■//—)■ 
III — )■ / and / — i- /// JJ J in terms of all probabilities p,. These 
transition rates must be equal, since DBC implies that transition rates for 
cyclic changes at equilibrium should not depend on their direction. 

Here for clarity we write the spins that are being changed in parentheses. 
Let us concentrate for a while on the transition (I) — t- (II). There are only 
two ways of choosing two consecutive spins that can affect a given spin: they 
are nearest neighbors either to the left or to the right, the probability of 
each choice is 1/L. In each case in order to obtain the desired change, one 
neighboring spin of the chosen pair must be changed while the other must 
not. For the pair ++ (written below in brackets) to the left from the spin 
being changed: 

... (+)[++](+) ... 

i (B.4) 

(+)[++](-) 

the resulting probability is pi (probability of changing spin to the right of 
the chosen pair) times p'^^ = 1 — pi (probability of not changing spin to the 
left of the chosen pair). Let us recall that pi's are defined by the formulae 

(IT])-(jl]) and p[ = 1 — Pi. Similarly, for the choice of the pair (again in 

brackets) to the right from the spin being changed: 

... (+) + +(+)[--]-... 

i (B.5) 

... (+) + +(-)[—]-... 

the resulting probability is p2 (probability of changing spin to the left of the 
chosen pair) times p'l = 1 ~ pi (probability of not changing spin to the right 
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of the chosen pair). Thus the overall probability of change (I) — )■ (II) in one 
step reads 

Wi^n = {pip'i+P2p'i)/L. (B.6) 
In this way we calculate the following transition probabilities between: 

Wn^iii = iP2P'i+PiP'2)/L, 

Wiii^i = pI/l, 

Wu^i = {P2P[+P1P[)/L, (B.7) 

Wni^n = {pip[+p2P2)/L, 

Wi^iii = pI/L. 



Let us now write DBC for transitions between any pair of the states 
considered and after taking their product we obtain: 

Wi^n qii Wii^iii qui Wiii^j qj = (B.8) 
= Wi^iii qui Win qii Wii^i qi, 

where qx denotes the probability that the system is at the state X at equi- 
librium. Hence, if qiqnqiii 7^ 0, 

Wi^ii Wn^iii Wiii^i = Wi^iii Wiii^ii Wii^i. (B.9) 



Inserting probabilities f lB.6p and f lB.7p into flB.9p and rearranging the terms, 



we obtain the condition: 

PliPl+ P2)iP2- Pl)[p'lipl+ P1P2+ pI) + P'2P1P2] =0. (B.IO) 

Taking into account that Pi > and p^ > 0, the above equation can be 
satisfied if and only if 

pi = 1 or pi = p2. (B.ll) 



By the same argument, for another choice: 



(/): 


••• + --(+) + +(+) 


- - + 


(//): 


••• + --(+) + +(-) 


- - + 


(///) : 


••• + --(-) + +(-) 


- - + 



(B.12) 
(B.13) 
(B.14) 



12 



we obtain conditions: 

P2 = 1 or pi = p2. (B.15) 

Taking the conjunction of the conditions (IB. lip and (IB. 150 we finally restrict 
the admissible choice of pi and p2 to 

Pi=P2- (B.16) 



In order to obtain conditions for p^ and p4 we consider transitions between 
the states: 



(///) 

and between the states: 

(/) 

(11) 
(III) 



+ -(+) + -(-) + - + ••• 
+ -(+) + -(+) + - + ... 

+ -(-) + -(+) + - + ••• 



••• + + -(+) + -(-) + 
... + + -(+) + -(+) + 

••• + + -(-) + -(+) + 



The resulting conditions are 

P3 = 1 or p3 = Pi 
for the first set of the states, and 

P4 = 1 or P3 = Pi 
for the second set, respectively. This leads to: 

P3 = Pi- 



(B.17) 
(B.18) 
(B.19) 



(B.20) 
(B.21) 
(B.22) 



(B.23) 



(B.24) 



(B.25) 
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